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Abstract. An integrable dispersionless KdV hierarchy with sources (dKdVHWS) 
is derived. Lax pair equations and bi-Hamiltonian formulation for dKdVHWS are 
formulated. Hodograph solution for the dispersionless KdV equation with sources 
(dKdVWS) is obtained via hodograph transformation. Furthermore, the disper- 
sionless Gelfand-Dickey hierarchy with sources (dGDHWS) is presented. 
PACS number: 02. 03. Ik 



1 Introduction 

In recent years, research in the dispersionless hierarchies has become quite active (see, for 
example, PQ-^I] and references therein). Dispersionless hierarchies arise as the quasiclas- 
sical limit of the original dispersionfull hierarchies The operators in the Lax equations 
for dispersionfull hierarchy are replaced by phase functions for dispersionless hierarchy, 
commutators are replaced by Poisson brackets and the role of Lax pair equations by disper- 
sionless Lax pair equations. The dispersionless hierarchies have Hamiltonian formulation 
[HE] and many other aspects [UOEE], and several methods of solutions of dispersionless 
hierarchies have been formulated [3 IHl El EE HI] • 

The soliton equations with self-consistent sources (SESCS) are another type of inte- 
grable models and have important physical applications p2]-j2H]- There are some ways 
to derive the SESCS, for example, the Mel'nikov way [121 CHI EH EH ED] and the Leon 
approach [13 EH HHJ . In recent years, SESCS were studied based on the constrained flows 
of soliton equations which are just the stationary equations of SESCSsjTB] HI]. There 
are several methods for solving the SESCSs, for example, the inverse scattering method 
[T£| EH ED] , the matrix theory [21] , the d method and gauge transformation [221 EH] , and 
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the Darboux transformation [TT1 124[ 12*51 126| . By treating the variable x as the evolution pa- 
rameter and t as the 'spatial' variable respectively, and introducing Jacobi-Ostrogradiski 
coordinates, the SESCS has a t-type Hamiltonian formulation [2*7| I28j. 

This paper is devoted to the integrable dispersionless KdV hierarchy with sources 
(dKdVHWS). Considering the asymptotic expansion of the wave function, and taking the 
dispersionless limit of the KdV hierarchy with sources (KdVHWS), we can deduce the 
dKdVHWS. The Lax pair equations of the dKdVHWS can be deduced by the dispersion- 
less limit of the Lax pair equations of the KdVHWS. Similar to the dKdV hierarchy, the 
dKdVHWS has bi-Hamiltonian formulation and can be solved via hodograph transfor- 
mation. The Gelfand-Dickey hierarchy with sources (GDHWS) is the integrable general- 
ization of the Gelfand-Dickey hierarchy, and the corresponding integrable dispersionless 
hierarchy, i.e. the dGDHWS, can be deduced. 

This paper is organized as follows: in Section 2 we review some definitions and re- 
sults about the KdVHWS. In Section 3 we derive the dKdVHWS as well as its Lax pair 
equations by taking the dispersionless limit of the KdVHWS. In Section 4, we construct 
the bi-Hamiltonian formulation of dKdVHWS. Then we derive the hodograph solution 
for dispersionless KdV equation with sources (dKdVWS) in Section 5. In Section 6, we 
deduce the integrable dispersionless Gelfand-Dickey hierarchy with sources (dGDHWS). 
Some conclusion is made in Section 7. 



2 The KdV hierarchy with sources 



We first briefly review some definitions and results about the KdV hierarchy with sources 
in the framework of Sato theory. Given a pseudo-differential operator (PDO) of the form 

m 

L = d 2 + u (2.1) 

where d = J^, u = u(x, t), t = (£ 3 , t 5 , . . . ), and the wave function ip = ip(x, t), consider 
the Lax pair 

bifi = (2.2a) 



= B 2m -i^, (2.2b) 

where A is a parameter, B 2m -i = (Lt°t^) + > m — L 2, 3, . . . , and (A)+ here standing for 
the differential part of A. The compatibility condition of (2.2a) and (2.2b) gives rise to 

dL 

tt: = [5 2m _i,L], (2.3) 
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which is the well-known KdV hierarchy [29]. As was shown in 29j, the KdV hierarchy 
could be written as bi-Hamiltonian systems 

u t2m _ 1 = B —^ = B 1 —^, m = l,2,..., (2.4) 



where J#2> m -i = f h^m-idx is a functional of u, B = d, B\ = + ud + \u x are 

8u 



Hamiltonian operators, and "? m ~ 1 is the Euler-Lagrange derivative of the Hamiltonian 



J#2m-i defined as 



/ h 2m -idx = 2_^{-dy 

J l — n 



k <-> lb 2m-l 

k=0 



(2.5) 



Let us consider l(i 

n 

Bzm-x = B 2m -i + } j ij}kd~ 1 <j>k, m = 1,2,3,..., (2.6) 
fc=i 

where vpk = i>k{ x )t) and = <pk(x,t) satisfying 

Lifj k = \ k ipk, L*4> k = \ k 4>k, k = 1, . . . , n, (2.7) 

here L* = (—d) 2 + u = L is the adjoint operator of L, and A& is a constant, = 1, . . . , n. 
Then the KdV hierarchy with sources (KdVHWS) PH EH 122 ED! can be defined as 

° L [B2m-i,L], (2.8a) 



dt2m-i 

Li/J k = Afc^fc, (2.8b) 

= Afc(/> fc (2.8c) 

= AV>, (2.9a) 

Vw.! = S 2m _ 1 V, (2.9b) 

namely, under (2.8b) and (2.8c), the compatibility condition of (2.9a) and (2.9b) gives 
rise to (2.8a). 



with the Lax pair given by 
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3 Dispersionless limit 



Following the procedure introduced in [21 El E3 [TO] , we could derive the dispersionless 
hierarchy by taking the dispersionless limit of the initial system. Taking T = et, X = ex, 
and thinking of u(x, t) = u(*> = U(X, T) + 0(e) as e — > 0, L in (2.1) changes into 

L e = e 2 d 2 x + u(j, ~) = e 2 d\ + U(X, T) + 0(e), (3.1) 

where d x — It can be proved [3] that 

C = a e (L e ) =p 2 + U (3.2) 

satisfies 

^T 2m _i = {#2m-l,£}, (3.3) 

where cx e denotes the principal symbol [2|, the bracket {, } is the Poison bracket defined 
in 2D 'phase space' (p, X) as 

, , _ OA dB OAdB . . 

and B%m-\ = (C~^~) + now refers to nonnegative powers of p. We define (3.3) as the 
dispersionless KdV (dKdV) hierarchy and the first few equations are expressed as 

U Tl = U x , (3.5a) 

Ut 3 = \UU X , (3.5b) 

U n = ^-U 2 U X , (3.5c) 
35 

U T7 = -U 3 U X ,... (3.5d) 
lb 

As was shown in [3], equation (3.5b) has tri-Hamiltonian formulation as 

m SH 5 _ SH 3 3_ 5 Hi , , 

r/ T3 = ^ = = jIV^, (3.6) 

where 

X>x = 2d Xl V 2 = Ud x + d x U, V 3 = U 2 d x + d x U 2 , 
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// / UdX, H 3 = ^ I U 2 dX, #5 = ~ / U :i <JX. 



For the general case, define 



= —^-Tr£^ , (3.7) 
2m — 1 



where TrA = J ResAdX, and ResA is the residue of the general Laurent polynomial of 
the form A = 2~2t^oo a i(X)p l } i.e. the coefficient of the p -1 term, then the Hamiltonians 
(3.7) are in involution with respect to any of the three Poisson brackets 

{H 2m ^, Hu-fr = J dX 6 -^^V l 6 -^ = 0, i = 1, 2, 3, (3.8) 

and dKdV hierarchy (3.3) have tri-Hamiltonian formulation 

U T2 ^vj*?* 1 =V 2 5 -^= (2m-l)(2m-3) ^ 2m _ 3 
T2m " 1 1 <J*7 2 5U (2m -2) 2 3 <^7 ' ' ' V ; 

It was also shown in [3] that the solution of (3.5b) can be described through the 
implicit form 

U = f(X + ~UT 3 ), (3.10) 

where / is an arbitrary function. 

In what follows, we derive the dispersionless KdV hierarchy with sources (dKdVHWS). 
By taking T = et, X = ex, (2.8) change into 

n XT XT 

eL eT 2m -i = [Be(2m-1) + Y] ^k (~ , ~) {^X ) "Vfc ( — ,—), L e] , (3.11a) 

f— ' e e e e 

k=l 

L € M-,-)=^kM-,-), (3-llb) 
e e e e 

L:0 fc (=|,|) = A fc fc (^,|), (3.11c) 



where 5 e ( 2m -i) = (L £ 2 ) + . 

Similar to the dispersionless KP case in (21 El ITU] , we consider the following WKB 
asymptotic expansion of i>k{—, 7) and </>&(—, 7), & = 1, 2, . . . , n, 

M*> 7) ~ exp{ g(I ' r ' 6 A = A ^ + (3 kl + 0(e)}, e - 0, (3.12a) 
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M*, 7) - exp{- S{X > T > e X - Xk) + (3 k2 + 0(e)}, e - 0, 



it can be calculated that 

M* ,7)WVfc(f .7) 
= e^^Kedx)- 1 + (A 5 (x,T,A = A fc ) + 0(e))(ed x )- 2 

+ (( J^(X, T, A = A,)) 2 + 0(e))(ed x r 3 + ...}, 

as e — > 0. Therefore, we have 



P-Pk 



where 

Vk = e /3 fel+ /3 fe2? pk = £_ S (X,T,\ = A fc ). 
Taking the principal symbol of both sides of (3.11a), we have 

and the dispersionless limit of (3.11b) and (3.11c) lead to 

d 

pl + U = A fc , -^(pkVk) = 0. 

Under (3.16), it can be found that 

{ — ,£>} = -2f fe ,x, 
P-Pk 

therefore, the dispersionless limit of (3.11), i.e. dKdVHWS, reads 

n 

u T 2m -i = {B 2m -i,£}-2^2v k ,x, m = 1,2,3,..., 

k=l 

p\ + U = A fc) 



ttv(r^) = 0. 



Integrating (3.18c) and taking A™ as the integral constants, then we can eliminate Vk in 
(3.18a) and rewrite dKdVHWS in another form 

V^_ l={ B^C)-^( 7 =M=) x (3.19) 

If we take the dispersionless limit of (2.9), we will obtain the Lax pair equations of 
dKdVHWS (3.18) as 

p 2 + U = A, (3.20a) 



n 

Pt^ = (B 2m -i + Yl -^-) > ( 3 - 20b ) 



namely, under (3.18b) and (3.18c), the compatibility condition of (3.20a) and (3.20b) gives 
rise to (3.18a). We can eliminate v\. and pk in (3.20b) and rewrite (3.20) in another form 

p 2 + U = A, (3.21a) 



\ 771 

pt^ = + n — rr k n — jw) > ( 3 - 21b ) 

v f^PV^k - U- (A fe - U)Jx 

which are the Lax pair equations of (3.19). 

We give two examples in the following, the first one is the dispersionless KdV equation 
with sources (dKdVWS) 

U T ^luU x -,t{-^=) x , (3,2) 
with Lax pair equations given by 

p 2 + u = A, (3.23a) 

^ ={f + lu^t ^*^^ . (3,3b) 
And the second example is the dispersionless KdV(5) equation with sources (dKdV(5)WS) 

, r n > 3 
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with Lax pair equations given by 

p 2 + U = A, 



(3.25a) 



m = (/ + + ^ + ± (3 , 5b ) 



4 Hamiltonian formulation of dKdVHWS 

It is well known that the KdV hierarchy has bi-Hamiltonian formulation [21], the dKdV 
hierarchy has tri- Hamiltonian formulation 3J, or even further, quasi-Hamiltonian for- 
mulation [1 , and the KdVHWS has bi-Hamiltonian formulation [27] . Motivated by the 
Hamiltonian formulation of the dKdV case [3], we would construct the bi-Hamiltonian 
formulation of the dKdVHWS (3.19). 

Let us firstly consider dKdVWS (3.22). Set 

A k = 2\ 2 k f V^k-U dX, B k = 2X k J VX k -U dX, (4.1) 

then by direct computation we have 

SA k nn / XI \ / XI 



VX k -U) \^X k -UJx 



<r\ /TT r, , a TT\f ^ k \ tfpx n ( Xi 

v i^T7 = ( ud x + d x U)[ = - = -2( 



su K A A J \ Vh=uJ (a* -co 3 / 2 y^^uJx 

Therefore, if we denote 

n „ 1 n 

H 3 = H 3 + J2Bk= / dX(-U 2 + 2j2Wh-U), 
H 5 = H 5 + Y,A k = / dX(-U 3 + 2j2^W*k-U), 

k=l J k=l 

then equation (3.22) can be written in two Hamiltonian forms 



(4.2a) 
(4.2b) 



For the dKdVHWS (3.19), denote 



n « 

H 2m -i + 2 K' 1 / VX k - U dX, (4.4) 



H2rn-1 

k=l 
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we can directly prove (see Appendix) that the Hamiltonians i?2m-i ; m = 1,2, . . . satisfy 
{H 2m ^, H 2 ^h = J dX 5 -^^V l 6 -^ = 0, i = 1, 2, (4.5) 
therefore, the dKdVHWS (3.19) have bi-Hamiltonian formulation 

^r 2ro -i = Pi = P 2 g , m = 1, 2, . . . . (4.6) 

5 Hodograph solution for dKdVWS 

In this section, using the hodograph transformation [7[ l9l ITU] . we will derive the hodograph 
solution for the the dKdVWS (3.22). Following [7] and let C% = B{U)U X , we can find 
from (3.22) that 

Making the hodograph transformation with the change of variables (X, T 3 ) — > ({7, T 3 ) and 
let X = X(£/, T 3 ), we have 

_ dX dU dX _ dX OX 

°-dT 3 -dum + dT 3 - du BUx + w*> (5 - 2) 

which implies that 

wr- B= ^ + £jA*- (5 - 3) 

It can be integrated as 



^2 

where F{U) is an arbitrary function of U, (5.4) gives an implicit solution of (3.22). 

When F = and Ai = A 2 = • • • = X n = 0, dKdVWS (3.22) degenerates to dKdV 
equation (3.5b), and (5.4) degenerates to the rational solution of (3.5b), U = — |^ . 

When F / and is convertible, then solution of dKdVWS (3.22) can be written 
through the implicit form 

/3 n A 2 \ 

U=F^X + (- 2 U-Y,^-^)T 3 ), (5.5) 

which is similar to (3.10). 
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6 The dispersionless Gelfand-Dickey hierarchy with sources 



The well-known Gelfand-Dickey hierarchy with sources (GDHWS) 16J is defined as 

dL n 

— = [B m ,L] = lB m + ^4>kd- 1 <Pk,L], (6.1a) 
0tm fc=i 

Lipk = \ki>k, (6.1b) 

L*(p k = h<Pk, (6.1c) 

where 

L = 8 N + u N ^ 2 8 N - 2 + • • • + u x d + u , (6.2) 

u = (mtv-2, • • • , u ) T , Ui = Ui(x, t), i = 0, 1, . . . , N - 2, t = (£ 2 , h, . . . ), B m = [(L^) m ]+, 
L* = (—d) N + (—8) n ~ 2 un-2 + • ■ ■ + {—d)u\ + uo is the adjoint operator of L, \ k is a 
constant, k = 1, . . . ,n, ifj k = ip k (x, t) and <p k = <fik( x , t), and the Lax pair is given by 

Lip = (6.3a) 

A m = Bm*Pi (6.3b) 

namely, under (6.1b) and (6.1c), the compatibility condition of (6.3a) and (6.3b) gives 
rise to (6.1a). 

Following the procedure given above, we can derive the dispersionless Gelfand-Dickey 
hierarchy with sources (dGDHWS). Taking T = et, X = ex, and let u k (—, — ) = U k (X, T) + 
0(e) as e — > 0, then (6.1) change into 

eL eTm = [B m + ^4 k (-,-)(ed x )- 1 (f> h (-,-),L e ] } (6.4a) 

L e i/j k (-, 7) = hMj, -), (6.4b) 
L:0 fc (^,^) = A fe fc (^,^), (6.4c) 



where 



(ed x ) N + -)(ed x ) N ~ 2 + ■■■ + Ul (-, -)ed x + u (-, -) 

e e e e e e 

(ed x f + (U X - 2 {X, T) + 0{e)){ed x ) N - 2 + ... 
+(U 1 (X,T) + 0(e))ed x + U (X,T) + 0(e), (6.5) 
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and B em = [(Li l ) m ] + . Consider the following WKB asymptotic expansion of ^( — , — 
and fc (f , ~), fc = 1, . . . ,n, 



^(-, -) ~ exp{ g(X ' r ' A ~ Afc) + (3 kl + 0(e)}, e - 0, 
e e e 



(6.6a) 



-) ~ ex P {- S{X ^ X - Xk) + ? k2 + 0(e)}, e -> 0, 



e e e 
then the principal symbol [2j of (6.4a) arises 



dC 



k=l 



P-Pk 



where 



£ = a e (L e ) = p N + t/ iV -2P 7V - 2 + • • • + U lP + U 



B m — ^{Bem), an d v fc — e /3fcl+/3fc2 , = -^S(X, T, X = X k ) are obtained by 

X 7\ 



P-Pk 



The dispersionless limit of (6.4b) and (6.4c) give rise to 



p% + u N - 2 pr z + ---+u lPk +u = X k , 



JV-2 



(6.6b) 



(6.7) 



(6.8) 



(6.9) 



(6.10a) 



d t dC 



dX \ Vk ' dp 



p=Pk 



0, 



(6.10b) 



(6.7) together with (6.10a) and (6.10b) give rise to the dispersionless Gelfand-Dickey 
hierarchy with sources (dGDHWS) 



dC 

dfl. 



(6.11a) 



p% + C/jv-2Pfc~ 2 + • • • + U lPk + U = X k , 



JV-2 



(6.11b) 



d 



dC 



dX \ Vh ' dp 



p=Pk 



0, 



whose Lax pair equations are given by 

p N + U N . 2 p N ~ 2 H h Uip + U = X, 



(6.11c) 



(6.12a) 
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PT m =(B m + J2^^) , (6.12b) 

namely, under (6.11b) and (6.11c), the compatibility condition of (6.12a) and (6.12b) gives 
rise to (6.11a). 

Under (6.11b) and (6.11c), it can be found by a tedious computation that 

{^^, C) = a N . 2 p N ~ 2 + a N ^p N ~ 3 + ■ ■ ■ + oq, (6.13) 
P-Pk 

where a N _ 2 = -Nv ktX , a N _ 3 = -N(v k p k ) x , and for % = 4, . . . , N, 

i-3 i-2 

a N -i = -v k Y,3P i k lu N + i-i + i,x - J2( N - l)M7 2 ~ l )xU N ^ - N(pl 2 v k ) x . (6.14) 
i=i 1=2 

When iV = 2, we have C = p 2 + U, and (6.13) is the same as (3.17). 

Similar to the dKdVHWS, the dGDHWS possess bi-Hamiltonian formation and their 
solutions can be obtained via hodograph transformation. 



7 Conclusion 



We derive dKdVHWS by taking the dispersionless limit of KdVHWS, meanwhile, Lax 
pair equations of the dKdVHWS can be obtained by taking the dispersionless limit of the 
corresponding dispersionfull equations. We have constructed the bi-Hamiltonian formula- 
tion of the dKdVHWS, and have obtained the implicit solutions of the dKdVWS via the 
hodograph transformation. For the generalization case, we have deduced the dGDHWS 
which also possess bi-Hamiltonian formulation and can be solved via the hodograph trans- 
formation. 
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Appendix 

Here we give the proof of involution relation of the Hamiltonians H 2m _i (4.4). 
Since 

{#2m-i, H 2 i-i}i = {H 2m -i, H 2 is}2 = — {H 2 i- 3 , H 2m _ 1 } 2 = —{H 2 i_ 3 , H 2m+ i}i 

= {H 2m+ i, H 2 i_ 3 } 1 = ■ ■ ■ = {H 2m+2 i_ 3 , Hi}i, 
it suffices to prove that for any m > 1, {H 2m -i, = 0. We can directly calculate that 

{H 2m -i, Hi}i 



5U 5U 



SB 



dX 



Since all SH $u~ 1 , rn = 1,2,... are of the form cll s , where c are constants, and s G Af, it 
suffices to prove 

j 7 ^TW^=VV^¥ ) dX = k,l = l,...,n. (U) 
For (J), let F l = t/X^U, then U = A, - F, 2 , = -2F,F, )JC and 



J Vk^u(k-u) J 



dF h 



which are all zero for the reason that ^'^p - are rational polynomials of Ff and so 
^ Xl ~3 ^ dFi are total differentials. 

For (II), when k — I, then (II) obviously holds; when k ^ I, let Gi = Xi — U, then 
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U x = -Gi >x , and 




set = tan 2 H kh then 1 + = sec 2 H kh d(^^) = 2tanH H sec 2 H kl dH kh and 
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